The Generalized Lyapunov Theorem and its Application to Quantum Channels 
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We give a simple and physically intuitive necessary and sufficient condition for a map acting on a 
compact metric space to be mixing (i.e. infinitely many applications of the map transfer any input 
into a fixed convergency point). This is a generalization of the "Lyapunov direct method". First we 
prove this theorem in topological spaces and for arbitrary continuous maps. Finally we apply our 
theorem to maps which are relevant in Open Quantum Systems and Quantum Information, namely 
Quantum Channels. In this context we also discuss the relations between mixing and ergodicity (i.e. 
the property that there exist only a single input state which is left invariant by a single application 
of the map) showing that the two are equivalent when the invariant point of the ergodic map is 
pure. 

I. INTRODUCTION 

Repetitive applications of the same transformation is the key ingredient of many controls techniques. In quantum 
information processing [ll they have been exploited to inhibit the decoherence of a system by frequently perturbing 
its dynamical evolution @, [3, 0, S @| {Bang-Bang control) or to improve the fidelity of quantum gates @ by means 
of frequent measurements {quantum Zeno-effect [8|). Recently analogous strategies have also been proposed in the 
context of state preparation [§, EfJ, EH, El, El, fl^.Tla| and quantum communication [lj| El, El, EH- In Refs. [ll], EH 
for instance, a homogenization protocol was presented which allows one to transform any input state of a qubit into 
a some pre-fixed targ et state by repetitively coupling it with an external bath. A similar thermalization protocol was 
discussed in Ref. [131 ] to study the efficiency of simulating classical equilibration processes on a quantum computer. In 
Refs. (l4l . [l5| repetitive interactions with an externally monitored environment were instead exploited to implement 
purification schemes which would allow one to extract pure state components from arbitrary mixed inputs. An 
application to quantum communication of similar strategies has been finally given in Refs. [m [l7l EH El] where 
sequences of repetitive operations were used to boost the efficiency of quantum information transmission along spin 
chains. 

The common trait of the proposals [1, El, EH, E3, EU EH, El, EH E3, El, El] is the requirement that repeated 
applications of a properly chosen quantum operation r converges to a fixed density matrix X* independently from the 
input state x of the system, i.e. 
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with "o" representing the composition of maps. Following the notation of Refs. [H, Hl| we call Eq. |1| the mixing 
property of r. It is related with another important property of maps, namely ergodicity (see Fig. [T|). The latter 
requires the existence of a unique input state xq which is left invariant under a single application of the map 1 , i.e., 
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Figure 1: Schematic examples of the orbits of a ergodic and a mixing map. 



Definition (f2j) may sound unusual for readers who are familiar with a definition of ergodicity from statistical mechanics, where a map 
is called ergodic if its invariant sets have measure or 1. The notion of ergodicity used in the case of a discrete time evolution of a 
quantum system is different l2Cj 1 12 otl . Here, the map r is not acting on a measurable space but on the compact convex set of quantum 
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Ergodicity and the mixing property are of high interest not only in the context of the above quantum information 
schemes. They also occur on a more fundamental level in statistical mechanics and open quantum systems [23L l24l] . 
where one would like to study irreversibility and relaxation to thermal equilibrium. 

In the case of quantum transformations one can show that mixing maps with convergency point x* are also ergodic 
with fixed point xo — x* ■ The opposite implication however is not generally true since there are examples of ergodic 
quantum maps which are not mixing (see the following). Sufficient conditions for mixing have been discussed both 
in the specific case of quantum channel [13l . |2& l2q ] and in the more abstract case of maps operating on topological 
spaces [22j]. In particular the Lyapunov direct method [22| allows one to prove that an ergodic map r is mixing if 
there exists a continuous functional S which, for all points but the fixed one, is strictly increasing under r. Here 
we strengthen this criterion by weakening the requirement on S: our generalized Lyapunov functions are requested 
only to have limiting values S(T n (x))\ n ^ oc which differ from S(x) for all x ^= xq. It turns out that the existence of 
such S is not just a sufficient condition but also a necessary condition for mixing. Exploiting this fact one can easily 
generalize a previous result on strictly contractive maps [20| by showing that maps which are asymptotic deformations 
(see Definition [6]) are mixing. This has, unlike contractivity, the advantage of being a property independent of the 
choice of metric (see however [2l| for methods of finding "tight" norms). In some cases, the generalized Lyapunov 
method permits also to derive an optimal mixing condition for quantum channels based on the quantum relative 
entropy. Finally a slightly modified version of our approach which employs multi- central Lyapunov functions yields 
a characterization of (not necessarily mixing) maps which in the limit of infinitely many applications move all points 
towards a proper subset (rather than a single point) of the input space. 

The introduction of a generalized Lyapunov method seems to be sound not only from a mathematical point of view, 
but also from a physical point of view. In effect, it often happens that the informations available on the dynamics of 
a system are only those related on its asymptotic behavior (e.g. its thermalization process), its finite time evolution 
being instead difficult to characterize. Since our method is explicitly constructed to exploit asymptotic features of 
the mapping, it provides a more effective way to probe the mixing property of the process. 

Presenting our results we will not restrict ourself to the case of quantum operations. Instead, following [22| we will 
derive them in the more general context of continuous maps operating on topological spaces [261 ] . This approach makes 
our results stronger by allowing us to invoke only those hypothesis which, to our knowledge, are strictly necessary 
for the derivation. It is important to stress however that, as a particular instance, all the Theorems and Lemmas 
presented in the paper hold for any linear, completely positive, trace preserving map (i.e. quantum channel) operating 
on a compact subset of normed vectors (i.e. the space of the density matrices of a finite dimensional quantum system). 
Therefore readers who are not familiar with topological spaces can simply interpret our derivations as if they were 
just obtained for quantum channels acting on a finite dimensional quantum system. 

The paper is organized as follows. In Sec. [TT] the generalized Lyapunov method along with some minor results are 
presented in the context of topological and metric spaces. Then quantum channels are analyzed in Sec. IIIII providing a 
comprehensive summary of the necessary and sufficient conditions for the mixing property of these maps. Conclusions 
and remarks end the paper in Sec. IIVI 

II. GENERALIZED LYAPUNOV THEOREM 
A. Topological spaces 

In this section we introduce the notation and derive our main result (the Generalized Lyapunov Theorem). The 
properties of Hausdorff, compact and sequentially compact topological spaces will be used [26|. For the sake of 
readability their definitions and their relations are given in the caption of Fig. O 

Definition 1. Let X he & topological space and let r : X — > X be a map. The sequence x n = r n (x), where r n is 
a short-hand notation for the n— fold composition of r, is called the orbit of x. An element € X is called a fixed 
point of r if and only if 

t(x*) = . (3) 



states. A perhaps more intuitive and equivalent definition of ergodicity based on the time average of observables is given by Lemma [5] 
of the Appendix. 
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r is called ergodic if and only if it has exactly one fixed point, r is called mixing if and only if there exists a convergence 
point x* G X such that any orbit converges to it, i.e. 

lim x n = Vx G X . (4) 

n — >oo 

Remark. Here we use the usual topological definition of convergence, i.e. limn^oo x n — x„ if and only if for each 
open neighborhood 0{x lr ) of x* only finitely many points of the sequence are not in 0(a;*). This clearly depends 
on the topology, and there may exist many different points to which a sequence converges to. For example, in the 
trivial topology of X where the only open sets are X and the empty set, any sequence is convergent to any point. 
On the other hand the uniqueness of the convergence point can be enforced by requiring the topological set X to be 
Hausdorff [26[ (see Fig. [2] for an explicit definition of this property). 

A direct connection between ergodicity and mixing can be established as follows. 
Lemma 1. Let r : X — > X be a continuous mixing map on a topological Hausdorff space X . Then r is ergodic. 




Figure 2: Relations between topological spaces [2(|. Hausdorff topological spaces have the property that any two distinct points 
of the space can be separated by open neighborhoods: for these sets any convergent sequence converges to a unique point 
of the set. Compact topological spaces are such that any open cover of the set has a finite sub-cover. Sequentially compact 
topological spaces are those for which the Bolzano- Weierstrass theorem holds, i.e. every sequence has a convergent subsequence. 
Any compact metric spaces is Hausdorff, compact, and sequentially compact. The space of density matrices on which quantum 
channels are defined, is a compact and convex subset of a normed vectors space (the space of linear operators of the system) 
which, in the above graphical representation fits within the set of compact metric spaces. 



Proof. Let ir* be the convergence point of r and let x € X arbitrary. Since r is continuous we can perform the limit 
in the argument of r, i.e. 

t(x*) = r( lim T n (x) ) = lim r" +1 (x) = x*, 

\n — >oo / n — >oo 

which shows that x* is a fixed point of r. To prove that it is unique assume by contradiction that r possesses a 
second fixed point y» ^ x». Then limn^oo T n (y*) = j/* ^ i t , so t could not be mixing (since the limit is unique in a 
Hausdorff space - see Fig. [2]). Hence r is ergodic. □ 

Remark. The converse is not true in general, i.e. not every ergodic map is mixing (not even in Hausdorff topological 
spaces). A simple counterexample is given by r : [—1, 1] — > [— f , 1] with t{x) = —x and the usual topology of R, which 
is ergodic with fixed point 0, but not mixing since for x ^ 0, r"(x) = (— l) n x is alternating between two points. A 
similar counterexample will be discussed in the quantum channel section (see Example [1} . 

A well known criterion for mixing is the existence of a Lyapunov function [22| . 

Definition 2. Let r : X — > X be a map on a topological space X. A continuous map S : X — > R is called a (strict) 
Lyapunov function for r around x* € X if and only if 

S (t(x)) > S(x) \/x ^ x*. 
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Remark. At this point is is neither assumed that ir* is a fixed point, nor that r is ergodic. Both follows from the 
theorem below. 

Theorem 1 (Lyapunov function). Let r : X — > X be a continuous map on a sequentially compact topological space 
X . Let S : X — > R be a Lyapunov function for r around x*. Then r is mixing with the fixed point cc*. 

The proof of this theorem is given in [12] • We will not reproduce it here, because we will provide a general theorem 
that includes this as a special case. In fact, we will show that the requirement of the strict monotonicity can be much 
weakened, which motivates the following definition. 

Definition 3. Let r : X — ► X be a map on a topological space X. A continuous map S : X — > R is called a generalized 
Lyapunov function for r around x* £ X if and only if the sequence S (r"(i)) is point-wise convergent for any x G A* 
and S fulfills 

SUrc) = lim S(r n (x)) ^ S(x) Vx ^ x». (5) 

n — >oo 

In general it may be difficult to prove the point-wise convergence. However if S is monotonic under the action of r 
and the space is compact, the situation becomes considerably simpler. This is summarized in the following Lemma. 

Lemma 2. Let r : X — > X be map on a compact topological space. A continuous map S : X — > R which fulfills 

S{t(x))^S(x) Vx£X, (6) 

and 

S*(x) = Km S(j n (x)) > S(x) Vx^x». (7) 

n — >oc 

for some fixed x* £ X is a generalized Lyapunov function for r around x* . 

Proof. It only remains to show the (point-wise) convergence of S (r™(x)). Since S is a continuous function on a compact 
space, it is bounded. By Eq. ^ the sequence is monotonic. Any bounded monotonic sequence converges. □ 

Corollary 1. Let t : X — > X be a map on a compact topological space. A continuous map S : X —> R which fulfills 

S(t(x))^S(x) Vx£X, 

and 

S (t n (x)) > S(x) Vx^x*, 
for some fixed N £ N and for some x» £ X is a generalized Lyapunov function for r around x* . 
Remark. This implies that a strict Lyapunov function is a generalized Lyapunov function (with N = 1). 
We can now state the main result of this section: 

Theorem 2 (Generalized Lyapunov function). Let r : X — > X be a continuous map on a sequentially compact 
topological space X. Let S : X — > R be a generalized Lyapunov function for r around x*. Then t is mixing with fixed 
point x* . 

Proof. Consider the orbit x n = t" (x) of a given x £ X . Because X is sequentially compact, the sequence x n has a 
convergent subsequence (see Fig. [5]), i.e. lim^oo x„ fc = x. Let us assume that i/i, and show that this leads to a 
contradiction. By Eq. (J5]) we know that there exists a finite N £ N such that 

S{t n {x))^S{x). (8) 

Since t n is continuous we can perform the limit in the argument, i.e. limfe^oo t n (x nk ) = t n (x). Likewise, by 
continuity of S we have 

lim S (x„ fc ) = S(x), (9) 

k — >oc 

lim S (x N+rik ) = lim S (r (x„J) = S(t n x), (10) 

k — >oo k — >oo 

where the second equality stems from the continuity of the map S and r . Because S is a generalized Lyapunov 
function, the sequence S (x n ) is convergent. Therefore the subsequences ([9]) and (fT0|) must have the same limit. We 
conclude that S(t n x) = S(x) which contradicts Eq. ([5]). Hence x = x*. Since we have shown that any convergent 
subsequence of r"(x) converges to the same limit x*, it follows by Lemma 3] of the Appendix that r n (x) is converging 
to x». Since that holds for arbitrary x, it follows that t is mixing. □ 



and on the other hand 
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There is an even more general way of denning Lyapunov functions which we state here for completeness. It requires 
the concept of the quotient topology . 

Definition 4. Let r : X — > X be a map on a topological space A". A continuous map S : X — > R is called a multi- 
central Lyapunov function for t around T C A" if and only if the sequence 5 (r"(x)) is point-wise convergent for any 
x £ X and if 5 and t fulfill the following three conditions: S is constant on J 7 , t(T) C JF, and 

S*(x) = lim 5 (r n (x)) ^ Six) Vx 4 T. 

n — >oo 

For these functions we cannot hope that the orbit is mixing. We can however show that the orbit is "converging" 
to the set T in the following sense: 

Theorem 3 (Multi-central Lyapunov function). Let t : X — > X be a continuous map on a sequentially compact 
topological space X. Let S : X — > M be a multi-central Lyapunov function for r around T . Let p : X — > XjT be 
the continuous mapping into the quotient space (i.e. <p(x) = [x] for x £ X\J- and <p(x) = [J 7 ] for x £ J 7 ). Then 
f : X jT — > X jT given by f([x]) = ip (t (p^ 1 ([%]))) is mixing with fixed point [J-]. 

Proof. First note that f is well defined because ip is invertible on X/T\[T] and t(T) C J 7 , so that = [J 7 ]. 

Since X is sequentially compact, the quotient space X jT is also sequentially compact. Note that for O open, 
f -1 (0) = ip (t -1 (p^ 1 (O))) is the image of <p of an open set in X and therefore (by definition of the quotient 
topology) open in X jT. Hence f is continuous. The function 5([a;]) : XjT — * XjT given by 5([x]) = S(p^ 1 ([x})) 
is continuous and easily seen to be a generalized Lyapunov function around [T]. By Theorem [5] it follows that f is 
mixing. □ 



B. Metric spaces 



We now show that for the particular class of compact topological sets which posses a metric, the existence of a 
generalized Lyapunov function is also a necessary condition for mixing. In this context the convergence of a sequence 
is defined with respect to the distance function d(-, ■) : X x X — > R on the space, so that for instance Eq. ((3]) requires 
linin^oo d(x n ,x») = 0. 

Theorem 4 (Lyapunov criterion). Let r : X — > X be a continuous map on a compact metric space X. Then t is 
mixing with fixed point x* if and only if a generalized Lyapunov function around x* exists. 

Proof. Firstly, in metric spaces compactness and sequential compactness are equivalent, so the requirements of The- 
orem [2] are met. Secondly, for any mixing map r with fixed point x*, a generalized Lyapunov function around x* is 
given by S(x) = d(x*, x). In fact, it is continuous because of the continuity of the metric and satisfies 

lim 5(r"(x)) = d(x*,x„) = ^ d(x*,x) = S(x), 

n — >oo 

where the equality holds if and only x = x*. We call d(x*,x) the trivial generalized Lyapunov function. □ 

Remark. In the above Theorem we have not used all the properties of the metric. In fact a continuous semi-metric 
(i.e. without the triangle inequality) would suffice. 

The trivial Lyapunov function requires knowledge of the fixed point of the map. There is another way of charac- 
terizing mixing maps as those which bring elements closer to each other (rather than closer to the fixed point). 

Definition 5. A map r : X — > X is on a metric space is called a non- expansive map if and only if 

d(r(x),T(y)) < d(x,y) Vx,yeX, 

a weak contraction if and only if 

d(r{x), r(y)) < d(x, y) Vx, y £ X, x ^ y, 

and a strict contraction if and only if there exists a k < 1 such that 

d(r(x),r(j/)) < kd(x,y) Vx,y £ X . 

Remark. The notation adopted here is slightly different from the definitions adopted by other Authors [13, HH 
who use contraction to indicate our non-expansive maps. Our choice is motivated by the need to clearly distinguish 
between non-expansive transformation and weak contractions. 
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We can generalize the above definition in the following way: 

Definition 6. A map r : X — > X on a metric space is called an asymptotic deformation if and only if the sequence 
d(r™ (x) , T n (y)) converges point-wise for all x,y £ X and 

lim d(r n (x),T n (y)) + d(x,y) Vx,y£X,x^ y. 

n — >oo 

Remark. Let r : X — > X be a non-expansive map on a metric space X, and let 

lim d(T n (x),T n (y)) < d(x,y) Vx, y £ X, x =/= y. 

71 — >00 

Then r is an asymptotic deformation. Any weak contraction is an asymptotic deformation. 

Theorem 5 (Asymptotic deformations). Let t : X — > X be a continuous map on a compact metric space X with at 
least one fixed point. Then r is mixing if and only if r is an asymptotic deformation. 

Proof. Firstly assume that r is an asymptotic deformation. Let be a fixed point and define S(x) = d(x*, x). 
lim S{r n {x)) = lim d(x*, T n (x)) 

n — >oo n — >oo 

= lim d(r n (x*), T n (x)) ^ d(x*, x) — S(x) \fx ^ x*, 

n — >oo 

hence S(x) is a generalized Lyapunov function. By Theorem [2] it follows that r is mixing. Secondly, if r is mixing, 
then 

lim d(r n (x),T n (y)) = d(x*,x*) = ^ d(x,y) Vx, y £ X, x ^ y, 

n — *oo 

so r is an asymptotic deformation. □ 

Remark. Note that the existence of a fixed point is assured if r is a weak contraction on a compact space 27], or if the 
metric space is convex compact [28] . As a special case it follows that any weak contraction r on a compact metric space 
is mixing. This result can be seen as an instance of Banach contraction principle on compact spaces. In the second 
part of the paper we will present a counterexample which shows that weak contractivity is only a sufficient criterion 
for mixing (see Example[2]). In the context of quantum channels an analogous criterion was suggested in [20l. I25I] which 
applied to strict contractions. We also note that for weak and strict contractions, the trivial generalized Lyapunov 
function (Theorem [4]) is a strict Lyapunov function. 



III. QUANTUM CHANNELS 

In this Section we discuss the mixing properties of quantum channels [l[ which account for the most general evolution 
a quantum system can undergo including measurements and coupling with external environments. In this context 
solving the mixing problem (fT]) is equivalent to determine if repetitive application of a certain physical transformation 
will drive any input state of the system (i.e. its density matrices) into a unique output configuration. The relationship 
between the different mixing criteria one can obtain in this case is summarized in Fig. [3] 

At a mathematical level quantum channels correspond to linear maps acting on the density operators p of the system 
and satisfying the requirement of bein g co mpletely positive and trace preserving (CPT). For a formal definition of 
these properties we refer the reader to |29l . [30. l3lll : here we remind only that a necessary and sufficient conditions to 
being CPT is to allow Kraus decomposition~"f30l] or. equivalently, Stinespring dilation [32 ] . Our results are applicable 
if the underlying Hilbert space is finite dimensional. In such regime there is no ambiguity in defining the convergence 
of a sequence since all operator norms are equivalent (i.e. given two norms one can construct an upper and a lower 
bound for the first one by properly scaling the second one) . Also the set of bounded operators and the set of operators 
of Hilbert-Schmidt class coincide. For the sake of definiteness, however, we will adopt the trace-norm which, given 
the linear operator 6 : TL — > H, is defined as ||6||i = Tr]V 9^9] with Tr[- • •] being the trace over TL and O^ being 
the adjoint of 9. This choice is in part motivated by the fact [33j that any quantum channel is non-expansive with 
respect to the metric induced 2 by || • \\i (the same property does not necessarily apply to other operator norms, e.g. 
the Hilbert-Schmidt norm, also when these are equivalent to || • ||i). 



2 This is just the trace distance d(p,a) = ||p — oj|i. 
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Figure 3: Relations between the different properties of a quantum channel. The red text indicates the new results obtained in 
this paper, and the black text indicates formerly known results which we reviewed. 



We start by showing that the mixing criteria discussed in the first half of the paper do apply to the case of quantum 
channels. Then we will analyze these maps by studying their linear extensions in the whole vector space formed by the 
linear operators of TL. Similar questions also arise in the context of finitely correlated states, where one investigates 
the decay of correlations in space (rather than in time) [34l ]. 



A. Mixing criteria for Quantum Channels 



Let TL be a finite dimensional Hilbert space and let S(TL) be the set of its density matrices p. The latter is a convex 
and compact subspace of the larger normed vector space C(H) composed by the linear operators 9 : TL — > TL of TL. 
From this and from the fact that CPT maps are continuous (indeed they are linear) it follows that for a quantum 
channel there always exists at least one density operator which is a fixed point [l3| . It also follows that all the results 
of the previous section apply to quantum channels. In particular Lemma [1] holds implying that any mixing quantum 
channel must be ergodic. The following example shows however that it is possible to have ergodic quantum channels 
which are not mixing. 

Example 1. Consider the qubit quantum channel r obtained by cascading a completely decoherent channel with a NOT 
gate. Explicitly r is defined by the transformations r(|0)(0|) = |1)(1|, = |0)(0|, and t(|0)(1|) = t(|1)(0|) = 

with |0), |1) being the computational basis of the qubit. This map is ergodic with fixed point given by the completely 
mixed state (|0)(0| + |l)(l|)/2. However it is trivially not mixing since, for instance, repetitive application of r on 
1 0> <0 1 will oscillate between |0)(0| and |1)(1|. 

Theorems [5] implies that a quantum channel r : S(TL) — > S(TL) is mixing if and only if it is an asymptotic deforma- 
tion. As already pointed out in the introduction, this property is metric independent (as opposed to contractivity) . 
Alternatively, if the fixed point of a quantum channel is known, then one may use the trivial generalized Lyapunov 
function (Theorem 0]) to check if it is mixing. However both criteria depend on the metric distance, which usually 
has no easy physical interpretation. A more useful choice of the is the quantum relative entropy, which is defined as 

H(p,a) = Trp(\ogp-loga). (11) 

As discussed in (3o| , the quantum relative entropy is continuous in finite dimension and can be used as a measure of 
"distance" (though it is not a metric). It is finite if the support of p is contained in the support of a. To ensure that 
it is a continuous function on a compact space, we consider the case when a is faithful: 

Theorem 6 (Relative entropy criterion). A quantum channel with faithful fixed point p* is mixing if and only if the 
quantum relative entropy with respect to p* is a generalized Lyapunov function. 

Proof. Because of Theorem [5] we only need to prove the second part of the thesis, i.e. that mixing channels admit 
the quantum relative entropy with respect to the fixed point, S(p) = H(p, p»), as generalized Lyapunov function. 
Firstly notice that the quantum relative entropy is monotonic under quantum channels [36| . Therefore the limit 
S*(p) = lrniji—xx, S (t t1 (jo)) does exist and satisfies the condition S*(p) ^ S(p). Suppose now there exists a p such that 
S*(p) — S(p). Because r is mixing and S is continuous we have 

S(p) = S*{p) = lim S(T n (p)) = S( P ,) = 0, 

71 — ► OO 

and hence H(p, p#) — 0. Since H(p, a) = if and only if p — a it follows that S is a Lyapunov function around p*. □ 
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Corollary 2 (Unital channels). A unital channel is mixing if and only if the von Neumann entropy is a generalized 
Lyapunov function. 

Another sufficient condition for mixing is weak contractivity. As already mentioned in the previous section, un- 
fortunately this not a necessary condition. Here we present an explicit counterexample based on a quantum channel 
introduced in Ref. [l3j |. 

Example 2. Consider a three-level quantum system characterized by the orthogonal vectors |0), |1), |2) and the quan- 
tum channel r defined by the transformations r(|2)(2|) = |1)(1|, r(|l)(l|) = r(|0)(0|) = |0)(0|, and r(|i)(j|) = for 
all i 7^ j . It's easy to verify that after just two iterations any input state p will be transformed into the vector |0)(0|. 
Therefore the map is mixing. On the other hand it is explicitly not a weak contraction with respect to the trace norm 
since, for instance, one has 

|| r(|2)(2|) - r(|0)(0|) ||r = |l |1)(1| - |0>(0| || x = || |2)(2| - |0)(0| ||i , 

where in the last identity we used the invariance of || • ||i with respect to unitary transformations. 



B. Beyond the density matrix operator space: spectral properties 

Exploiting linearity quantum channels can be extended beyond the space S(H) of density operators to become 
maps defined on the full vector space L(H) of the linear operators of the system, in which basic linear algebra results 
hold. This allows one to simplify the analysis even though the mixing property ([I]) is still defined with respect to the 
density operators of the system. 

Mixing conditions for quantum channels can be obtained by considering the structure of their eigenvectors in the 
extended space C{7i). For example, it is easily shown that the spectral radius (37j of any quantum channel is equal 
to unity |l3|, so its eigenvalues are contained in the unit circle. The eigenvalues A on the unit circle (i.e. |A| = 1) are 
referred to as peripheral eigenvalues. Also, as already mentioned, since S(H) is compact and convex, CPT maps have 
always at least one fixed point which is a density matrix (l3j . A well-known connection between the mixing properties 
and the spectrum is given by the 

Theorem 7 (Spectral gap criterion). A quantum channel is mixing if and only if its only peripheral eigenvalue is 1 
and this eigenvalue is simple. 

Proof. The "if" direction can be found in linear algebra textbooks (see for example (37l . Lemma 8.2.7]. Now let us 
assume that r is a mixing quantum channel with fixed point p». Let be a generic operator in C{TL). Then can 
be decomposed in a finite set of non-orthogonal density operators 3 , i.e. 6 = J2e c ePe> with pi g S(H) and a complex. 
Since Tr [pi] = 1, we have have Tr [O] = ^2 e eg. Moreover since r is mixing we have linin^oo T n (pi) = p* for all £, 
with convergence with respect to the trace-norm. Because of linearity this implies 

lim r" (9) = Vqp»= Tr [6] p* . (12) 

I 

If there existed any other eigenvector 0* of r with eigenvalue on the unit circle, then linin^oo t™(0») would not 
satisfy Eq. ([12]). □ 

The speed of convergence can also be estimated by [l]| 

\\r n {p)-pA\i ^C N n N n n , (13) 

where N is the dimensionality of the underlying Hilbert space, k is the modulus of the second largest eigenvalue of r, 
and Cat is some constant depending only on N and on the chosen norm. Hence, for n» N the convergence becomes 
exponentially fast. As mentioned in [2fl|, the criterion of Theorem [7] is in general difficult to check. This is because 
one has to find all eigenvalues of the quantum channel, which is hard especially in the high dimensional case. Also, 
if one only wants to check if a particular channel is mixing or not, then the amount of information obtained is much 
higher than the required amount. 



3 To show that this is possible, consider an arbitrary operator basis of C(7i). If N is the finite dimension of fi. the basis will contain 
N 2 elements. Each element of the basis can then be decomposed into two Hermitian operators, which themselves can be written as 
linear combinations of at most N projectors. Therefore there exists a generating set of at most 27V 3 positive operators, which can be 
normalized such that they are quantum states. There even exists a basis (i.e. a minimal generating set), but in general it can not be 
orthogonalized. 
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Example 3. As an application consider the non mixing CPT map of Example^ One can verify that apart from the 
eigenvalue 1 associated with its fixed point (i.e. the completely mixed state), it possess another peripheral eigenvalue. 
This is X = — 1 which is associated with the Pauli operator |0)(0| — |1)(1|- 

Corollary 3. The convergence speed of any mixing quantum channel is exponentially fast for sufficiently high values 
ofn. 

Proof. From Theorem [7] mixing channels have only one peripheral eigenvalue and it is simple. Therefore the derivation 
of Ref. [Hi applies and Eq. CCH holds. □ 

This result should be compared with the case of strictly contractive quantum channels whose convergence were 
shown to be exponentially fast along to whole trajectory [20l[25|. 



C. Ergodic channels with pure fixed points 



An interesting class of ergodic quantum channel is formed by those CPT maps whose fixed point is avure density 
matrix. Among them we find for instance the maps employed in the communication protocols of Refs. [la. ITtI ll8L ll9T| 
or those of the purification schemes of Refs. [3, EH- We now show that within this particular class, ergodicity and 
mixing are indeed equivalent properties. 

Theorem 8 (Purely ergodic maps). Let IV'iXV'il be the pure fixed point of an ergodic quantum channel t. It follows 
that t is mixing. 

Proof. We will use Theorem [7] showing that 1S the only eigenvector of r with peripheral eigenvalue. Assume 

in fact that S C(H) is a eigenvector of r with peripheral eigenvalue, i.e. 

r(e) = e 4V e. (14) 

From Lemma|6]of the Appendix we know that the density matrix p = V 00t / g, with g — Tr V 0t0 > 0, must be a 

fixed point of r. Since this is an ergodic map we must have p = \ipi)(ipi\. This implies = g\if}i)(ip2\, with \1jj2) some 
normalized vector of Tt. Replacing it into Eq. (fT4|) and dividing both terms by g yields t (1P2D = e IC, ° | -i/^i ) (^2 1 and 

\{i>M^ 1 }^ 2 \)\^ 2 )\ = i. 

Introducing a Kraus set {K n } n of r and employing Cauchy-Schwartz inequality one can then write 
1 = |^i|r(|^)^ 2 |)|V 2 )| = |^i|if„|Vi>^ 2 |^2>| 



where we used the fact that \tpi) is the fixed point of r. Since r is CPT the quantity (falTdfa) (faDlfa) is upper 
bounded by 1. Therefore in the above expression the inequality must be replaced by an identity, i.e. 

Since r is ergodic, we must have 1^2) (^2 1 = IV'iXV'iI- Therefore oc which shows that IV'iXV'il is the only 

peripheral eigenvalue of r. □ 

An application of the previous Theorem is obtained as follows. 

Lemma 3. Let Mab — Ma <8 Is + 1a ® Mb be an observable of the composite system H.a *S> T~Ib and r the CPT 
linear map on Ha of Stinespring form [32| 

T{p) = Tx B [U{p®\ct>) B {<t>\)U^ , (15) 

(here Trx [■ • •] is the partial trace over the system X, and U is a unitary operator of H.a ® 'Hb)- Assume that 
[Mab,U] — and that \4>) B is the eigenvector corresponding to a non- degenerate maximal or minimal eigenvalue of 
Mb- Then t is mixing if and only if U has one and only one eigenstate that factorizes as \v) a ® \4>)b- 
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Proof. Let p be an arbitrary fixed point of r (since r is CPT it has always at least one), i.e. Ttb \U (p® |</>)s(0|) U'] = 
p. Since Mab = Ma + M B is conserved, and Tr^ [Map] = Tr^ [M^r(p)], the expectation value of M B is unchanged. 
Hence system B must remain in the state with maximal/minimal eigenvalue, which we have assumed to be unique 
and pure, i.e. 

u( P ® \<f>)B{<f>\)tf = P®\<t>) B {<i>\ =► [u, P ®\<t>) B {<i>\] =o. 

Thus there exists a orthonormal basis {|i«fc)} fc of 7^ <S>TLb diagonalizing simultaneously both U and p® |</))s(0|. We 
express the latter in this basis, i.e. p® \<P)b{4\ — ^2kPk\ u k)(uk\, and perform the partial trace over subsystem A to 
get 

\<P)b(4>\ = ^PfcTr^ [\u k )(u k \] . 

k 

Hence Tr^ [|itfc)(ufc|] = \4>)b(4>\ f° r au k, and \u k ) must be factorizing, 

K) = \vk)A ® \<P)b- (16) 

If the factorizing eigenstate of U is unique, it follows that p = \v){v\ for some \v) and that r is ergodic. By Theorem 
[5] it then follows that r is also mixing. If on the other hand there exists more than one factorizing eigenstate, then all 
states of the form of Eq. (|16[) correspond to a fixed point p k = \v k ){y k \ and r is neither ergodic nor mixing. □ 

The case discussed in Lemma[3]is a generalization of the CPT map discussed in Ref. [l6[ in the context of spin-chain 
communication. There TLa and TLb represented two distinct part of a chain of spins coupled through Heisenberg-like 
interactions: the latter including the spins controlled by the receiver of the message, while the former accounting 
for all the remaining spins. Assuming the system to be initially in the ground state (i.e. all spin down), the sender 
(located at one of the extremes of the chain) encodes her/his quantum messages (i.e. qubits) into superpositions 
of spins excitations which will start propagating toward to receiver (located at the other extreme of the chain) . In 
Ref. [IH it was shown that, by repetitively swapping the spins which are under her/his control with some ancillary 
spins prepared in the ground state, the receiver will be able to recover the transmitted messages. The key ingredient 
of such result is the fact that by applying the swapping operations the receiver is indeed removing all the excitations 
(and therefore the corresponding encoded quantum information) out that part of the chain which is not directly 
accessible to him/her (i.e. the part represented by TLa)- In its simplest version, the resulting transformation on TLa 
can be described by Eq. (TT5)) with U and Mab representing, respectively, the free evolution of the spins among two 
consecutive swaps and the z-component of the magnetization of the chain. Lemma [3] can then be used to provide an 
alternative proof of convergence of the protocol [1 61 ] showing that indeed repetitive applications of t will drive TLa 
toward a unique convergence point (i.e. the state with no excitation). 

IV. CONCLUSION 

In reviewing some known results on the mixing property of continuous maps, we derived a stronger version of the 
direct Lyapunov method. For compact metric spaces (including quantum channels operating over density matrices) it 
provides a necessary and sufficient condition for mixing. Moreover it allows us to prove that asymptotic deformations 
with at least one fixed point must be mixing. 

In the specific context of quantum channels we employed the generalized Lyapunov method to analyze the mixing 
properties. Here we also analyzed different mixing criteria. In particular we have shown that an ergodic quantum 
channels with a pure fixed point is also mixing. 
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APPENDIX 



Here we derive some Lemmas which are not correlated with each other but which are relevant in our discussion. 
Lemma S] discusses a property of sequentially compact topological spaces. Lemma [5] states a well known theorem (20j 
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which, in the context of normed vector spaces, shows the equivalence between the definition of ergodicity of Eq. (|4|) 
and its definition using time averages. Finally Lemma [5] discusses a useful property of quantum channels (sec also 

M)- 

Lemma 4. Let x n be a sequence in a sequentially compact topological space X such that any convergent subsequence 
converges to x*. Then the sequence converges to x*. 

Proof. We prove by contradiction: assume that the sequence does not converge to x*. Then there exists an open 
neighborhood 0{x*) of x* such that for all k there is a n k such that x nk ^ 0(x*). Thus the subsequence x Hk is 
in the closed space X\0(x*), which is again sequentially compact. x nk has a convergent subsequence with a limit in 
X\0(x*), in particular this limit is not equal to x». □ 

Lemma 5. Let X be a convex compact subset of a normed vector space, and let r : X — > X be a continuous map. If 
t is ergodic with fixed point x* , then 

1 - 

l im y T e {x)=x, . (17) 

n^oc n + 1 z — ' 

Proof. Define the sequence A n = ^j-j- Y^t=o T ^( x )- Let then M be the upper bound for the norm of vectors in X, i.e. 
M = sup^g^. ||x|| < oo. which exists because X is compact. The sequence A n has a convergent subsequence A nk with 
limit A. Since r is continuous one has lim^oo r(A nh ) = t(A). On the other hand, we have 

1 „- fc +i^ „„ / \h nk+ \x)\\ + IWI / 2M 



-(A nk ) - A n J = — \\T n « +1 (x) - x\\ < ^ LJJ1_JLJ1. <; 



n fe + 1 ' n fc + 1 n k + 1 ' 

so the two sequences must have the same limit, i.e. t(A) = A. Since r is ergodic, we have A = x« and lim,^^ A n — x* 
by Lemma |U □ 

Remark. Note that if r has a second fixed point y» 7^ 2;* , then for all n one has 53"=o T ^ (f * ) = ^* ' so ^1- 03 
would not apply. 

Lemma 6. Let r be a quantum channel and be an eigenvector of t with peripheral eigenvalue A = e lip . Then, given 
g = Tr V 0t0 > 0, i/ie density matrices p = V 00t/<7 and er = V Q^Q/g are fixed points of t. 

Proof. Use the left polar decomposition to write = g pU where U is a unitary operator. The operator pU is clearly 
an eigenvector of r with eigenvalue e l<fi , i.e. 

r(pC/) = A pU . (18) 

Hence introducing a Kraus set {-fC„},i of r [30j and the spectral decomposition of the density matrix p — y 1 . Pj\ipj){ipj \ 
with pj > being its positive eigenvalues, one gets 



A = Tr[r(pC/)[/t] = £ p^K^^UK^^t) 



where the trace has been performed with respect to an orthonormal basis {|<^)}^ of Ti. Taking the absolute values of 
both terms gives 

|A| = ^PMKn^mUKlU^t)] 



= y/Tr[T(p)]y/Tr[f(p)] = 1, 

where the inequality follows from the Cauchy-Schwartz inequality. The last identity instead is a consequence of the 
fact that the transformation f(p) = Ut(U' pU)U' is CPT and thus trace preserving. Since |A| = 1 it follows that the 
inequality (|19[) must be replaced by an identity. This happens if and only if there exist e such that 

y/Pj{(4n\K n \<l>j)Y = VPji^Klh) = e^^p-^UKlU^) , 
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for all j,£ and n. Since the \<pi) form a basis of Tt, and pj > this implies 



, ii) 



-a) 



for all n and for all the not null eigenvectors \4>j) of p. This yields 



r(p)C/ 



which, replaced in (| 18[) gives e 1 r(p) = e l¥> /?, whose only solution is 



e i</> xhcrefore r(/?) = p and p is a fixed 



point of t. The proof for a goes along similar lines: simply consider the right polar decomposition of O instead of the 
left polar decomposition. □ 

Corollary 4. Let r be an ergodic quantum channel. It follows that its eigenvectors associated with peripheral eigen- 
values are normal operators. 

Proof. Let G be an eigenoperator with peripheral eigenvalue e lip such that r (O) = e lip 8. By Lemma[6]we know that, 
given g — Tr V 8^8 the density matrices p = VO&/g and a = V&^d/g must be fixed points of r. Since the map 
is ergodic we must have p = <r, i.e. 88^ = 0*0. □ 
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